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Abstract 

We theoretically show that in a crystal with a helical lattice structure, orbital and spin magne¬ 
tizations along a helical axis are induced by an electric current along the helical axis. We propose 
a simple tight-binding model for calculations, and the results can be generalized to any helical 
crystals. The induced magnetizations are opposite for right-handed and left-handed helices. The 
current-induced spin magnetization along the helical axis comes from a radial spin texture on the 
Fermi surface. This is in sharp contrast to Rashba systems where the induced spin magnetization 
is perpendicular to the applied current. 


* murakami@stat.phys.titech.ac.jp 


1 




■ novel physics due to an interplay between electricity and magnetism 
1-4], skyrmion jiii-sj], and current-induced magnetization reversal [s]- 
12], makes it possible to control a magnetization with electric held. However, there is 


Recent discovery o 
such as multiferroics 


a well-known classical example; an electric current flowing through a solenoid induces a 
magnetization. In this paper, we theoretically propose a condensed-matter analogue of 
a solenoid. We consider a three-dimensional crystal with a helical lattice structure. For 
crystals with helical crystal structure such as Se or Te, one can define handedness similar to 
a solenoid. Keeping these crystals in mind, we propose a simple tight-binding model with 
handedness, which can grasp the essence of the physics of helical crystals. Then we show that 
orbital and spin magnetization along a helical axis is induced by an electric current along 
the helical axis. The induced magnetizations are opposite for right-handed and left-handed 
helices. In contrast to Rashba system, the spin texture on the Fermi surface is radial, and 
the current-induced spin magnetization along the helical axis in helical crystals is attributed 
to this radial spin texture. Our results can be generalized to any crystals without mirror 
and inversion symmetries and they would pave the way to spiutronics application of helical 
crystals. 

First, we consider an orbital magnetization in a helical crystal. Here, we introduce a 
three-dimensional tight-binding model with a right- or left-handed helical structures. The 
lattice structure of this model is composed of an infinite stack of honeycomb lattice layers 
with one orbital per site. We consider the honeycomb lattice as shown in Fig©, with 
bi = ax, bo = a/2(—x + \/3y), and b ,3 = —bi — b 2 , where a is a constant. The layers 
are stacked along the 2 -direction with the primitive lattice vector a 3 = cz, where c is the 
interlayer spacing. The Hamiltonian is 


H = ti c\cj + t 2 c\cj + A Cj, (1) 

(b> [bi * 

where H, t 2 and A are real and we set t\ > 0 for simplicity. The first term is a nearest- 
neighbor hopping term within the xy plane. The second term represents “helical” hoppings 
between the same sublattice in the neighboring layers. This term is different between a 
right-handed and a left-handed helices as shown in Fig©) and Fig©. In the right-handed 
helix (Fig©)), the direction of hoppings between A sites are ±(b, + cz), and that between B 
sites are ±(—b, : + cz). Similarly in the left-handed helix (Fig©) the direction of hoppings 
between A sites are ±(—b* + cz), and that between B sites are ±(bj + cz). This term breaks 
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inversion and mirror symmetries. The third term is a staggered on-site potential, where is 
+1 and —1 for the sites in A and B sublattices, respectively. The space group of our model 
is -P321, with C‘ 2 y and C :iz symmetries. Only when A = 0, the space group becomes /hi22, 
having additional and Cq z symmetries. 

We note that the Hamiltonian for the right-handed (left-handed) helix becomes the Hal¬ 
dane model on a honeycomb lattice [hij] by replacing +k z c (— k z c) with a Aharonov-Bohm 
phase (j) in the second-neighbor hopping in the Haldane model. We also note that the time- 
reversal symmetry is present in our model, while it is absent in the Haldane model for 
(j) 0. Hence the band structure of our model is easily obtained from that of the Haldane 


model 


13]. The Brillouin zone (BZ) and band structure of our model are shown in Figs. [2h 


and, b and c, respectively. In the BZ, the K (K 7 ) point is defined by k • bj = —2n/3 (2 tt/3) 
on the k z = 0 plane, and the H (H 7 ) point is similarly defined on the k z = 7r/c plane. 

Using this model, we calculate the orbital magnetization induced by an electric held 
along the helical axis. In the limit of zero temperature T —>• 0, it is calculated from the 


formula 


14 
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M orb — 2 


' BZ 


/„ so/nk m»tk T , (&F £rik)^?ik 
(27t) j L n 


( 2 ) 


where the integral is performed over the BZ, n denotes the band index, f n k is the distribution 
function for the eigenenergy £ n k, £p is the Fermi energy, and the prefactor 2 comes from the 
spin degeneracy. The orbital magnetic moment of the Bloch electrons is defined by 
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m„k = — i(e/2h) (delink \ x[H k — £„k] | <9k« n k) and the Berry curvature is defined by f2„k = 
i (delink | x |<9knnk), where \u n k) is the periodic part of the Bloch function. Because our 



FIG. 1. Helical lattice structure of the present model, a, One layer of the model, forming 
a honeycomb lattice. Dashed arrows denote vectors bi,b2, and b3. b, Hopping texture in the 
right-handed helix, c, Hopping texture in the left-handed helix. Red (blue) balls are hoppings 
between A (B) sites. 
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FIG. 2. Current-induced orbital magnetization, a, Brillouin zone of our model with high- 
symmetry points and that on k z = 0. b, c, Energy bands for the Hamiltonian ([!]) with = fi/3, 
and b A = H; c A = 3t\. The energy bands are indicated within 0 < k z < 7r/c. d, e, Orbital 
magnetization for several values of A in d the right-handed helix and e the left-handed helix. We 
set the parameters £2 = ii/3. 


model lacks inversion symmetry, Both m n k and f2 n k are allowed to have nonzero values for 
any k. However, in equilibrium, due to the time-reversal symmetry M orb is zero because in 
Eq. ([2]) the contribution from k and that from — k cancel each other. To induce the orbital 
magnetization, we apply an electric held along the z axis. For a metal it induces a charge 
current accompanied by nonequilibrium electron distribution, and orbital magnetization is 
expected to arise. Within the Boltzmann approximation, the applied electric held E z changes 
fn k into 


ink = fnk + eE z TV r 


df_ 

de 


£ ^nk 


(3) 


where /° k = f(e n k) is the Fermi distribution function in equilibrium, r is the relaxation 
time assumed to be constant, and v njZ = (l/h)de n k/dk z is the velocity in the z direction. 
Substituting f n k into Eq. ©, we obtain the current-induced orbital magnetization M or b, 
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which is along the z-axis by symmetry. Because of time-reversal symmetry, the /° k term in 
Eq. (J3J) does not contribute to M or b. In the T —> 0 limit, df/de has a sharp peak at Ef- 
Therefore, for a band insulator, the induced orbital magnetization M or b is zero. This is in 
contrast to multiferroics. Furthermore, the Berry-curvature term (second term in Eq. ((2])) 
always vanishes at T —> 0. 

Figure [2] d, e shows numerical results of Eq. (|2]) for several values of A. We set the 
parameters as to = ti/3. The band structure is easily obtained from the similarlity to the 


Haldane model 


13]. For |A| > 2ti, the bands have a finite gap. For y/3ti < |A| < 2 ti, the 


two bands overlap and the system is in the semimetal phase. For |A| < y/3t\, the two bands 
cross each other only on K-H line and K’-H’ line at A = — y/3t\ sin k z c and A = y/3ti sin k z c, 
forming Dirac cones at energies — (A 2 /3 tf), respectively. As is expected, when the 

Fermi energy lies in the energy gap (A = 3H, — 1\ < Ef < H in Figj2]d. e), the orbital 
magnetization is zero for T —^ 0. For metals (A = 0 and A = t\) 1 the current-induced 
orbital magnetization appears. In this case, the induced magnetization is largely enhanced 
around Ep — t\. This is attributed to an enhanced orbital magnetic moment m n k near 
the Dirac points appearing on the K-H or K’-H’ lines. The orbital magnetization in the 
left-handed helix (Figj2^) is exactly opposite to that in the right-handed one (Figj2]d). This 
dependence on handedness is similar to the solenoid, where an electric current generates a 
magnetic field. Thus, based on the semiclassical theory, the orbital magnetization induced by 
the current is attributed to a helical motion of a wavepacket. Here, we note that although 
most of the magnetoelectric effect involve the spin-orbit interaction, the current-induced 
orbital magnetization predicted here appears even without the spin-orbit interaction. 

Next, we consider current-induced spin polarization. To this end, we introduce the spin- 
orbit interaction into a tight-binding model with the helical structure: 


Hso — h T: c\cj + (iV3\/a) Y, c \{ sX d X ij + s v d^ j )c j 
(ij) (ij) 

+(iKJa)J24(s*d,° j +s‘d’’ j )c j 

[hi 

+i}K/c) c l sZ dijCj + A Y & C K ( 4 ) 

[hi * 


The first term is a spin-independent nearest-neighbor hopping term within the xy plane 
and the fifth term is a staggered on-site potential term. The other three terms represent 
spin-orbit interactions (s are the Pauli matrices in spin space). The second term is a spin- 


5 




dependent nearest-neighbor hopping term within the xy plane. The third and forth term 
involve spin-dependent helical hoppings between the neighboring layers as shown in FigsHJ) 
and c, where is the /c-component of a vector pointing from site j to site i. The space 
group of this model is the same as that of the model ([!]). 

Figure [3^ shows the band structure of the Hamiltonian (J4J) along some high symmetry 
lines. Figure |3jo shows the spin texture projected onto the xy plane, on the Fermi surface 
around the H point at £p = 0.68ti. In fact, the spin around the H point has not only 
the x?/-component but also the z-comp orient. Because of the spin-orbit interaction, the 
two spin-split Fermi surfaces appear, having the opposite spin orientations. Remarkably, 
unlike Rashba systems, the spin is oriented radially, and rotates once around the H point. 
Figure [3jc shows the spin texture between the K and H points at £p = 0.68fi. In this case, 
one of the spin-split Fermi surface is open. Nevertheless, the inner Fermi surface have a 
radial spin texture. This spin texture results from crystal symmetries. Namely, the K-H 
lines are three-fold rotation axis and hence the spin on these lines are parallel to these K-H 
lines. Furthermore, the absence of mirror symmetries is crucial for the radial spin textures; 
if a mirror plane including the xy plane or the z axis were present, the spins should be 
perpendicular to the mirror plane, and a radial spin texture would not appear. For example, 


Te and Se have helica’ 


a radial spin texture 


crystal structure, and as expected, they have been predicted to show 


23]. 


We now calculate current-induced spin magnetization in the present system. It is known 
that in Rashba systems the current-induced magnetization is perpendicular to the current, 


because of the tangential spin texture in the spin-split. Fermi surfaces 


24 


27}. Similarly, a 


spin magnetization is induced by an electric current in our model, whereas the magnetization 
is parallel to the current when the current is along the helical axis by symmetry, thanks to 
the radial spin texture. Within the Boltzmann approximation as explained in the Methods 
section, the z-component of M sp i n for several parameters is numerically calculated as shown 
in Figure [3]d. We fix X xy = 0.05H, = 0.05H, and A = 0.4H- As is expected, for a 

metal, the spin magnetization parallel to the current is induced by the current. Similar to 
the orbital magnetization, the spin magnetization in the left-handed helix and that in the 
right-handed helix are opposite. 

We have shown that in systems with helical structure, the orbital and spin magnetizations 
are induced by an electric current along the helical axis, using a simple tight-binding models. 
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FIG. 3. Current-induced spin magnetization, a, Energy band for the Hamiltonian with 
A = — 0.06fi, X xy = X z = 0.05ti, and A = 0.4ti. Dashed line indicates ef = 0.68ti. b, c, Fermi 
surface and spin texture. The parameters are the same as a. Arrows represent the spin projected 
onto each plane. The spin on the inner Fermi surface in b is drawn with three times larger scale, 
d, Spin magnetization for several values of A. The parameters are \ xy = 0.05ti, A^ = 0.05fi, and 
A = 0.4fi. We also set the temperature T = 0.03ti/A;B- 


When an electric field is applied along the helical axis, the orbital magnetization is induced 
by the orbital magnetic moment on the Fermi surface, as a consequence of broken inversion 
symmetry. Similarly to the Berry curvature, the orbital magnetic moment is enhanced near 
band crossings. Therefore, when the Fermi energy lies near the band crossing, the orbital 
magnetization is enhanced as well. We have also shown that, in the helical crystal with 
spin-orbit interaction, a current induces the spin magnetization along the helical axis due 
to the radial spin texture. The absence of mirror symmetry allows to have the radial spin 
texture, which is completely different from Rashba systems with a tangential spin texture. 

We have presented a toy model for a helical structure. From a symmetry viewpoint, 
the present effect of longitudinal current-induced magnetization appears for chiral crystals 
without inversion and mirror symmetries. Real helical crystals, such as Se and Te, also lack 
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inversion and mirror symmetries. From symmetry consideration, we expect that orbital and 
spin magnetizations along the helical axis can be induced by an electric current in these 
crystals by doping carriers. Magnitudes of the induced magnetizations can be estimated 
from our results. The induced orbital magnetization scales with e 2 E z rti/h 2 , which is about 
150 gauss for E z = 10 4 V/m, r = 10~ 12 s and t\ = 3eV. In Fig. [21 the maximum value is 
0.1 times the above scale, and this factor could be enhanced for more “helical” crystalline 
structure. On the other hand, the induced spin magnetization scales wtih eF^r/ie/ha 2 , 
which is about 7 gauss for E z = 10 4 V/m, r = 10 _12 s and a = 0.5nm. It is multiplied by 
a numerical factor in Fig. [21 approximately given by the ratio X/ti. Thus, the size of the 
spin-orbit coupling limits the size of the current-induced spin magnetization. On the other 
hand, the orbital magnetization does not require spin-orbit coupling, and it can be enhanced 
by appropriate choice of materials with helical crystal structure. Te and Se consist of weakly 
coupled helices, and therefore they may have large current-induced orbital magnetization. 

Our results also provide a new building block of spintronics and will pave the way to 
spintronics application of helical crystals. For example, consider a ferromagnet on a helical 
crystal. By injecting current into the helical crystal, the induced magnetization exerts a 
torque on the ferromagnet, thus leading to current-induced magnetization reversal. 


Methods: Details of the calculation of current-induced spin magnetization- Within 
the Boltzmann approximation, the electric field E z induces the spin magnetization as 


M spin = - 


X ( Unk 


eE z ry B 

h 

df/ k 


E 


r d 3 k df_ 

Ibz (2vt) 3 de 


^ ^nk 


dk z 


(Unk | 7 <8> S | Unk) > 


(5) 


where is the Bohr magneton and the electron spin ^-factor g ~ 2. In the present model dH), for 
A = 0 the spin magnetization always vanish, because of a cancellation between the contribution 
from (k x , k y , k z ) and that from (—k x , —k y , i r/c — k z ). This cancellation can be avoided by adding 
either spin-dependent hoppings within the xy plane (for example, the A term in Eq. ([2|)) or spin- 
independent hoppings between the neighboring layers (for example, the term in Eq. (fT[) ). 
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I. PROPERTIES OF THE HAMILTONIAN ( 1 ) IN K SPACE 


The Hamiltonian (1) for the right-handed helix can be written as 


ifk = 2 t 2 cos (k z c) cos(k • bj 


+ 

+ ti 


1 + cos(k • a x ) + cos(k ■ a 2 ) 
sin(k ■ a x ) + sin(k ■ a 2 ) cr v 


a 


A 2t 2 sm(/ig c) sin(k • bj) 


a 


(SI) 


where / is the 2x2 identity matrix and a 1 are the Pauli matrices acting on the sublattice 
degree of freedom. Vectors a x and a 2 , defined as a x = a/ 2(x+\/3y) and a 2 = a/ 2(— x+\/3y), 
are the primitive lattice vectors within the xy plane. As we noted in the main text, the 
Hamiltonian H k becomes the Haldane model on a honeycomb lattice [1] by replacing k~c 
with a Aharonov-Bohm phase 0 in the second-neighbor hopping in the Haldane model. 
Therefore, the band structure of the Hamiltonian (ISID is easily obtained from that of the 
Haldane model [1] by replacing 0 with k z c. By analogy with the Haldane model, we assume 
|t 2 /fi| < 1/3. For |A/t 2 | > 6, the bands have a finite gap. For 3^ < |A/£ 2 | < 6, the 
two bands overlap and are not degenerate; namely systems are in semimetal phases. For 
| A/i 2 | < 3 a/ 3, the two bands cross each other and form Dirac cones, only on K-H and K’-H’ 
lines in the BZ. On the K-H line, the two bands cross when A = — 3\^3t 2 sin k z c, while on 
the K’-H’ line the two bands cross when A = 3\/3 1 2 sin k z c, as shown in Fig. [ST] 


II. PROPERTIES OF THE ORBITAL MAGNETIC MOMENT 


^he orbital magnetic moment for the electronic state at k of the nth band is defined 


by j2|, 3] 


(f^k^)ik | X [-^k £nk] | ^k^nk) • (S2) 

n 

The properties of the orbital magnetic moment are similar to those of the Berry curvature |4J. 
If the system has time-reversal symmetry, the orbital magnetic moment satisfies 

m n (-k) = —m n (k). (S3) 

If the system has spatial inversion symmetry, the orbital magnetic moment satisfies 

m n (—k) = m n (k). (S4) 
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FIG. SI. Positions of the Dirac points as a function of A. The Dirac points (Weyl nodes) 
appear along the K-H line (solid curve) and the K’-H’ line (broken curve). Only the values of k z 
are shown. 


Therefore, unless a system has both time-reversal and inversion symmetries, m„(k) is in 
general nonzero. 

The orbital magnetic moment can be also written as a summation over the eigenstates. 
The ^-component of the orbital magnetic moment is calculated as 


ze 


m n Jk) = -- 

= Jim 

n 


du 


nk 


dk r 


du 


nk 


dk~ 


[Hk <^nk] 
[-^k ^nk] 


du 


nk 


dk„ 


du. 


nk 


du 


nk 


dk. 


R Im S 

m^n 


^ n k 


dH± 

dk x 


^mk 


y 

^mk 


9H k 

dk v 


dk„ 


^nk 


[H\t ^nk] 


du, 


nk 


dkrr 


(^mk ^nk) 


(S5) 


The final expression does not include derivatives of the eigenstates, which is useful for 
numerical calculations. Equation (lS5l) shows that m n (k) diverges when e n k = £ m k at certain 
values of k. Therefore, m n (k) is enhanced around degeneracy points. 

We now consider the Hamiltonian (IS 11) . The orbital magnetic moments of the two bands 
are identical as seen from Eq. (1S5I) . Figure IS2h shows the energy band of the Hamiltonian 
(ISII) on the K’-H’ line and Fig. IS2b shows the ^-component of the orbital magnetic moment 
on the K’-H’ line. Dirac points appear at k z c = sin _1 (l/-\/3) and n — sin" - 1 (1/a/3). We can 
see that the orbital magnetic moment is enhanced around the Dirac points and diverges at 


3 


























a 



b 


100 



-50 

-100 






r a 



K’ sirr 1 (14/3) TT-sin- 1 (1A/5) H' 


k z c 


FIG. S2. a, Energy band for the Hamiltonian (ISlj) along the K’-H’ lines. The parameters are 
i 2 = t±/3, and A = t\. b, The z-component of the orbital magnetic moment on the K’-H’ line. 

the Dirac points. The enhanced orbital magnetic moment around the Dirac points makes a 
large contribution to the orbital magnetization. 


III. K-SPACE REPRESENTATION OF THE HAMILTONIAN (4) 


The Hamiltonian H so for the right-handed helix can be written in k-space as 


H SO} k 


The 2x2 matrices IIa■ H Bl and H AB are 


H a H ab 
H i B Hb 


(S6) 


H a = A/ + A xy [2 sin(k • b 3 + k z c ) — sin(k • b 2 + k z c ) — sin(k • b 3 + k z c))s x 
+ V3A X y [sin(k ■ b 2 + k z c) — sin(k • b 3 + k z c)\s y 

+ 2A,[y: sin(k • b^ + k z c)\s z , (S7) 

i 

H b = —A/ — X xy [2 sin(—k ■ bi + k z c ) — sin(—k • b 2 + k z c ) — sin(—k • b 3 + k z c)\s x 
— \/3A :ry [sin(—k • b 2 + k z c) — sin(—k • b 3 + k z c)\s y 

+ 2A z [y^sin(—k • b ;: + k z c)\s z , (S8) 


H ab = H[1 + e _ik ' ai + e-‘ ka2 ]J 

fo 1 

+ ^-iA[e.- 4k ai - e~ ika2 ]s x - -i\[2 - e 


(S9) 
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where / is the 2x2 identity matrix and s l are the Pauli matrices acting on spin space. 
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